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Abstract 

> ' 

" We consider the microscopic solutions of the Boltzmann-Enskog equation discov- 

ered by Bogolyubov. The fact that the time-irreversible kinetic equation has time- 
reversible microscopic solutions is rather surprising. We analyze this paradox and 
show that the reversibility or irreversibility property of the Boltzmann-Enskog equa- 
tion depends on the considered class of solutions. If the considered solutions have 
£NJ . the form of sums of delta-functions, then the equation is reversible. If the considered 

solutions belong to the class of continuously differentiable functions, then the equa- 
tion is irreversible. Also, we construct the so called approximate microscopic solutions. 
These solutions are continuously differentiable and they are reversible on bounded time 
?-h ' intervals. 



This analysis suggests a way to reconcile the time-irreversible kinetic equations 
with the time-reversible particle dynamics. Usually one tries to derive the kinetic 
equations from the particle dynamics. On the contrary, we postulate the Boltzmann- 
Enskog equation or another kinetic equation and treat their microscopic solutions as 
the particle dynamics. So, instead of the derivation of the kinetic equations from the 
microdynamics we suggest a kind of derivation of the microdynamics from the kinetic 
equations. 

1 Introduction 

In the known works of Bogolyubov [1, 2J, the Boltzmann-Enskog kinetic equation is shown to 
have microscopic solutions. This result is surprising for two reasons. The first reason is that 
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the Boltzmann-Enskog equation is regarded as an approximation to the particle dynamics 
in the case of a dilute gas. This result shows that the Boltzmann-Enskog equation has exact 
microscopic solutions as well. Also the Vlasov kinetic equation is known to have microscopic 
solutions [3]. 

But, in contrast to the Vlasov equation, the Boltzmann-Enskog equation is time- 
irreversible and describes the entropy production. However, the microscopic dynamics is 
time-reversible. So, this result tells that an irreversible equation has reversible solutions. 
The aim of this report is to analyze this paradox. 

We show that the reversibility or irreversibility property of the Boltzmann-Enskog equa- 
tion depends on the considered class of solutions. If the considered solutions have the form 
of sums of delta-functions, then the equation is reversible. If the considered solutions belong 
to the class of continuous functions, then the equation is irreversible. 

Also, we construct the so called approximate microscopic solutions of the Boltzmann- 
Enskog equation. These solutions are continuous and they are reversible on bounded time 
intervals. The approximate microscopic solutions turn to the (exact) microscopic solutions 
in the limit. 

This result can be compared with the Kozlov's diffusion theorems for a collisionless contin- 
uous medium in a rectangular box [I] . The dynamics is described by the Liouville equation. 
If we consider a solution in the form of a delta-function, then the solution is periodic or 
almost periodic (this case is reduced to the single-particle dynamics). But if we consider an 
integrable solution, the spatial density weakly converges to the uniform distribution, as time 
indefinitely increases. So, the asymptotic dynamical properties of solutions depend on the 
considered class of functions. 

But anyway, the Liouville equation is time-reversible. Here we show that in the case 
of the Boltzmann-Enskog equation, not only the asymptotic properties of solutions, but 
also the reversibility or irreversibility property of them depend on the considered class of 
functions. 

Further, this analysis suggests a way to reconcile irreversible kinetic equations with the 
reversible microscopic dynamics. Kinetic equations describe the gas as a whole, while equa- 
tions of microdynamics describe the same gas as a collection of particles (atoms or molecules). 
If we believe that the properties of the g whole can be reduced to the properties of 

its particles (the position of reductionism), then we shall postulate the equations of micro- 
dynamics and try to derive the corresponding kinetic equations from them. 

An elegant derivation of the Boltzmann and Boltzmann-Enskog equations from the Li- 
ouville equation was proposed by Bogolyubov [5]. He uses the BBGKY (Bogolyubov-Born- 
Green-Kirkwood-Yvon) hierarchy of equations, the thermodynamic limit and some addi- 
tional assumptions (which do not follow from the Liouville equation or its initial data). 
This derivation leads to the divergences in high order corrections to the Boltzmann and 
Boltzmann-Enskog equations [Bj. 

Another derivation of the Boltzmann equation was proposed by Lanford [7] (see also 
[HIE])- He also uses the BBGKY hierarchy of equations, the Boltzmann-Grad limit, all as- 
sumptions are included in the initial conditions for the Liouville equation. But the derivation 
can be applied only to small times. On the contrary, the Boltzmann equation is interesting 
from the viewpoint of the large time asymptotics (the relaxation to the Maxwell distribution) 
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Dm- 

So, we have not a satisfactory derivation of the Boltzmann and Boltzmann-Enskog equa- 
tions from the equations of microdynamics. 

In this report we propose an alternative way to reconcile kinetic equations with equations 
of microdynamics. We postulate not the equations of microdynamics, but the kinetic equa- 
tions and treat their microscopic solutions as the particle dynamics. The establishment of 
the existence of microscopic solutions can be regarded as a kind of derivation of the particle 
dynamics from a kinetic equation. So, instead of the derivation of kinetic equations from 
the microdynamics we suggest a derivation of the microdynamics from kinetic equations. 

In this interpretation the "fundamental" laws of motion are not the laws of motion for 
individual particles (reductionism), but the laws of motion for the single-particle density 
function, which describes the system as a whole (holism). Laws of motion for the individual 
particles are approximations. 

The Boltzmann-Enskog equation describes the dynamics of a hard sphere gas. In the 
final section of the paper we consider other kinetic equations. Namely, we consider the 
Boltzmann-Enskog equation with an additional "Vlasov integral" (we call it the Boltzmann- 
Enskog-Vlasov equation) and the Boltzmann equation. The Boltzmann-Enskog-Vlasov 
equation describes the dynamics of hard spheres that additionally interact with some long- 
ranged interaction potential. This equation also has microscopic solutions and, hence, one 
can speak about a derivation of the Newton equation (or Hamilton equations) of micrody- 
namics from it. 

The Boltzmann equation has not microscopic solutions, but it can be regarded as an 
approximation of the Boltzmann-Enskog (or a Boltzmann-Enskog-like) equation in the case 
of a large number of particles and a negligible interaction radius of them (the Boltzmann- 
Grad limit). In this sense one can speak about a reconciliation of the microdynamics with 
the Boltzmann equation as well. 

Recently Igor Volovich proposed functional mechanics [TH [12] (see also [13], HU HH [T6] ) 
for solving the irreversibility problem. The functional mechanics states that the fundamental 
notion of mechanics is not a material point, but an (integrable) probability density function. 
Hence, the fundamental equation even for a single particle is not the Newton equation (or 
Hamilton equations), but the Liouville equation (or even the Fokker-Planck-Kolmogorov 
equation). A solution of the Newton equation can be regarded as a solution of the Liouville 
equation in the form of a delta-function (another interpretation: the Hamilton equations 
are the equations of characteristics for the Liouville equation). But delta-functions are not 
allowed to be considered. A particle with definite position and momentum is an approxima- 
tion, which is suitable if the integrable density function has a specific (delta-like) form. 

As we have mentioned above, the spatial density for an integrable solution of the Liouville 
equation in a box weekly converges to the uniform distribution. So, if we postulate not the 
Newton, but the Liouville equation for a single particle, we obtain a kind of relaxation to the 
equilibrium state for it. Thus, the paradox between the recurrence of the microdynamics and 
the convergence to the equilibrium state of the macrodynamics is eliminated: both dynamics 
exhibit the convergence to the equilibrium. 

We use the same idea: we notice that the solution of the Newton equation for the system 
of hard spheres can be regarded as a special type solution (in the form of the sum of delta- 
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functions) of the Boltzmann-Enskog equation. Then we postulate the Boltzmann-Enskog 
equation. The existence of microscopic solutions allows us to treat the particle dynamics as 
an approximation, which is suitable if the density function has a specific (delta-like) form. 

Also Igor Volovich obtained a kinetic equation for a finite number of particles [17]. In fact, 
the Boltzmann-Enskog equation is also a kinetic equation for a finite number of particles. 

I am very happy to dedicate this report to the 65th anniversary of Professor Igor Volovich, 
and I would like to wish him good health and further scientific successes. 



2 Boltzmann-Enskog equation and its microscopic so- 
lutions 

Let us consider a gas of hard spheres with diameter a > in some 3-dimensional domain. In 
order to avoid difficulties with taking a boundary of the domain into account, let us assume 
that the domain is the 3-dimensional torus T 3 . The dynamics in a rectangular box is known 
to be reduced to the dynamics on a torus. The Boltzmann-Enskog kinetic equation for this 
gas takes the form 

dt ~~ Vl W x {) 

where 

St/ = na 2 / (v 2 i,a)[f(r 1 ,v[,t)f(r 1 + aa,v' 2 ,t)-f(r 1 ,v 1 ,t)f(r 1 -aa,Vi,t)]dadv 2 . (2) 

J(v2i,a)>0 

Here f(rx,Vx,t) is a density of particles near the point r% G T 3 that have velocities near 
V\ G M 3 at the moment of time t. Let / be normalized on the volume V of the torus T 3 : 

f(n,vx,t)drxdvx = V. 

T 3 xR 3 

The expression St fx is called the collision integral, n > is a constant (the mean 
concentration of particles), Vx and v 2 are the velocities of the colliding particles just before 
the collision, v[ and v' 2 are the velocities of the particles just after the collision: 

v'x = Vx + <t(v2x,o-), ^ 
v' 2 = v 2 - (t(v 21 ,(t). 

Here v 2 x — v 2 — v\ , the vector a is a unit vector from the center of the second particle to the 
center of the first particle (so, a belongs to the 2-dimensional unit sphere: o G S 2 ), (-, •) is 
a scalar product. 

The Boltzmann-Enskog equation describes the irreversible dynamics of a hard sphere gas 
and the entropy production, a rigorous .ff-theorem for this equation (and for the more general 
Enskog-type equations) can be found in [T8l [T9] . The if-function for the Boltzmann-Enskog 
equation slightly differs from the .ff-function for the Boltzmann equation. 
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Now let us consider a function / which has the form of a sum of delta-functions: 

N 

f(n, Vl , t; r) = n l 6 ^ - r )^i - Wi(t, r)), (4) 

8=1 

where gi(t, T), tOi(£, T), . . . , qx(t, T), Wx(t, Y) are the positions and velocities of iV hard 
spheres at the moment t. Here n = y. The time dependence is defined by the initial 
positions and velocities 

T = (q 1 ,wl...,q N ,w° N ) 

and the laws of motion for hard spheres (i.e., law ([3]) for pairwise collisions and the free motion 
between the collisions; we neglect the collisions of three and more spheres simultaneously, 
because the corresponding initial phase points T have zero measure as well as initial phase 
points which lead to the infinite number of collisions in a finite time [20J). 

Proposition 1 (Bogolyubov). The (generalized) function f given by Op satisfies the 
Boltzmann-Enskog equation (TJp. 

A solution of form (j4]) is called a microscopic solution of a kinetic equation. This result 
is surprising for two reasons. The first reason is that the Boltzmann-Enskog is regarded as 
an approximation to the iV-particle dynamics in the case of a dilute gas (n is small). This 
result shows that this equation has exact microscopic solutions as well. The Vlasov kinetic 
equation is also known to have microscopic solutions [3]. 

But, in contrast to the Vlasov equation, the Boltzmann-Enskog equation is time- 
irreversible and describes the entropy production. This proposition tells that it has solutions 
/, which are defined by the time-reversible dynamics of iV hard spheres. In other words, both 
f(ri,Vi,t;T) and f(ri, v±, t; T) = f(ri,—vx,—t;T) are solutions of the Boltzmann-Enskog 
equation, despite of the fact that this equation is time-irreversible! 

Indeed, f(r 1 ,-v 1 ,-t;T) = f(n,vi,t; T), where T = (q^, -w^, . . . , q%, -w° N ). So, the 
function / also has form (TJJ and, hence, is a solution of ([T|). 

3 Rigorous formulation of the microscopic solutions of 
the Boltzmann-Enskog equation 

Let us explore the relation between the irreversible Boltzmann-Enskog equation ([1]) and 
its reversible microscopic solutions 01]). In this section we define solutions / of form (jlj) 
rigorously. 

Denote G = T 3 x M. 3 (the single-particle phase space), Q is the set of points 
(ri, vi, . . . ,rN, Vn) € G N such that |r*, — Tj\ > a for all i ^ j. Obviously, red, since 
hard sphere centers cannot be closer than a to each other. Let D(T) be a continuously 
differentiable function on Q and rapidly decrease in Vi, . . . , Vpj with all its first-order partial 
derivatives in r and v. Let also the following condition be satisfied: 

D(. ..,n,Vi,...,Ti - aa,vj, . . .) = D{. . . ,r;,t> 4 ', . ..,n - aa,v' p . . .) (5) 
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for all i 7^ j, n, Vi, Vj, a and other arguments (denoted as ". . .") such that (vj — Vi, a) > 0. 
Here (r^i^) and (r* — aa,Vj) are zth and jth pairs of arguments of D, the other arguments 
are the same in both sides of the equalities, v[ and v'j are defined by (J3]) (with indexes 1 
and 2 replaced by i and j). Note that we can regard D as a function on G N with D = if 
|rj — 7j-| < a for some i ^ j. 
Put by definition 

N 

(II " r )), = ^^(xi, . . . , x N ) (6) 

i=i 

where Xi = (r^fj), ^(t) = (qi(t),Wi(t)). Here 5 t is the dynamic flow, i.e., 
S i t N \x l ,...,x N ) = (Xi(t),...,X N (t)), where X;(t) = (Ri(t),Vi(t)) 6 T 3 x K 3 are the posi- 
tions and velocities of N hard spheres at the moment t if their positions and momenta at 
the moment t — were (xi, . . . , xn). We use the notation 

StD(xi, x N ) = D(St(xi, x N )). 

Here sf 5 can be regarded as an operator which maps the function D to the function with 
the shifted arguments. 

The meaning of definition ([6]) is the following. It is natural to define the expression 
(n.i=i ${ x i ~ T)), D(r)) as the value of D in a point T such that yi(t,T) = x^ for all 
i — 1, . . . , N. This is exactly S_pD(xi, . . . , Xn). 

If D has the meaning of the initial probability density function, then S_P D = D t has 
the meaning of the probability density function at the time t. 

Further, 

(5(xi - yi(t, T)), D(T)) = / SlPd(x 1 , . . . , x N ) dx l . . . dx^idx i+ i . . . dx N . 

Jg n ~ 1 

So, we define the left-hand side of the expression as the integral from D(T) over the initial 
values T such that yi{t, Y) = x^. This definition can be obtained by the formal integration 
of (JHD over all variables Xi, . . . ,Xn, except X {. Analogously, 

(*(*, - Vl (t. r))6(xj - vA t. r)), D) = f sifDi Xl . . . dx ,„ ldx , +1 . . . dXi „ ldX]+1 

Jg n ~ 2 

Denote 

(f(x 1 ,t;T),D(T)) = F 1 (x 1 ,t), 

{f{x 1 , t- T)f{x 2 , t; T), D(T)) = F 2 {x u x 2 ,t). 

F\ and F 2 are the single-particle and two-particle density functions correspondingly. If 
D(x\, . . . , xjv) is symmetric relative to the permutations of Xi (the particles are indistin- 
guishable), then 

Fi(x\,t) = V / Sip D(xi, . . . ,xn) dx 2 . . . dxN, 

Ja "- 1 (8) 



(7) 



F 2 (x u x 2 ,t) = V(l-h [ S%>D( Xl , 



, xtv) dx$ . . . dx^j. 
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Condition (jSJ) can be rewritten in terms of the two-particle density function F 2 : 



F 2 (n,v 1 ,r 1 - aa,v 2 ,t) = F 2 (ri,v[,r 1 - aa,v' 2 ,t) 



(9) 



for all 7*1, Vi, v 2 , o and t such that (t> 2 i, cr) > 0. 

Note that the operator is discontinuous, since the velocities change jump-like at the 
moments of collisions. But if condition (jSJ) is satisfied, then S_^._ D = S^+ D. 

It would be desirable to define / as a generalized function on some space of test functions. 
We can define / as a functional on T>(Q\dQ) (infinitely differentiable functions with compact 
support in Q\dfl where dQ is a boundary of Q). Equality (jSJ) is trivially satisfied: both sides 
of it are equal to zero, since the points with |rj — rj\ = a for some i 7^ j belong to dQ. But 
the conditions imposed upon a function from T>(Q\dQ) (infinite differentiability and equality 
to zero on the boundary) are too restrictive, we want / to be defined on a wider class of 
functions. 

We can define / as a functional on the space of continuously differentiable functions on 
fl that rapidly decrease in Vi, . . . , v^. Denote this space as SC l {VL). A topology on SC 1 (fl) 
can be specified by the following seminorms: 



D\ 



sup 

n 



'1 + v' 



3D 



dri 



sup(l + v 

n 



3D 



dvi 



1,2, 



N}, 



+ . . . + v 2 Nl and K = 1, 2, But in this case, we must impose 



Pk(D) = max{sup(l + v A 
n 

where D e SC\Q), v 2 = 
additional functional condition flH]). Let us regard / as a functional (generalized function) 
on the subspace of SC 1 ^) defined by additional condition Denote this subspace as 
CC\tt). 

The meaning of condition (j5J) (or (Q) is easy to understand if D is a probability density 
function. This is a continuity condition: since the collision of hard spheres is instantaneous, 
the two-particle density just after the collision should be equal to the density just before 
the collision. If this condition is not satisfied, then the function D t is, generally speaking, 
discontinuous in x\, . . . ,xn- So, we demand the continuity not only for the initial density 
function, but for the density function at every moment of time as well. Also, this condition 
means the continuity of D t in t. 

In fact, we can drop condition (JHJ). But this case is more sophisticated to analyze. For 
simplicity we will assume condition flH]). 

Now we are able to formulate proposition 1 as a rigorous assertion. 

Proposition 2. The (generalized) function f given by (fjp satisfies the Boltzmann-Enskog 
equation (Tjp. The both sides of the equation are understood as functionals on CC 1 (f2). 

Proof. The single-particle and two-particle density functions for the system of iV hard spheres 
are known to satisfy the following equation (the first equation of the BBGKY hierarchy for 
hard spheres, see [21]; condition (E]) is used in the derivation of this equation): 



dF 1 (r 1 ,v l ,t) 
dt 



-«i Vna / {v 2 i,o-)[F 2 {ri,v 1 ,r l + aa,v 2 ,t) 

- F 2 (ri,v 1 ,r 1 - ao-,v 1 ,t)]dadv 2 , (10) 
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where n = y. Here the argument t\ ± aa means n ± aa ± 0, since .^(ri, ^1, r 2, ^2, = as 
|ri — r 2 | < a. Due to ([7]) this equation can be written as 

r r dfjruvutir) + 9/(r l Pl> t;r) _ ^ } ^ = 
J ^ at ori J 

where St / is given by (jSJ). Since D(T) is an arbitrary function from CC 1 (f2), this equality 
means that / satisfies the Boltzmann-Enskog equation in the stated sense. □ 

It should be noted that the microscopic solutions of the Boltzmann-Enskog equation are 
understood slightly differently than the microscopic solutions of the Vlasov equation. In the 
case of the Vlasov equation, a generalized function in form is understood as a functional 
on T>(G) (infinitely differentiable functions on G with compact support). So, if A G T>(G), 
then 

(8(r - qi (t, F),v- w { (t, T)),A(r, v)) = A(q t (t, T), Wi (t, T)). 

Thus, test functions depend on r and v. But it is difficult to give a sense to the expressions 
f(r\,vi, t)f(ri — aa, V2,t) and -£ in the Boltzmann-Enskog equation as functionals on T>(G). 
For this reason, they are defined on CC l (Q). In particular, test functions depend on T e Q. 
In this case, all terms can be well-defined. 

It could be seemed that the microscopic solutions of the Boltzmann-Enskog equation are 
just a paraphrase of equation flTDJ). However, the following theorem reveal another view on 
the microscopic solutions. 

4 Approximate microscopic solutions 

Theorem 1. Consider the function 

N 

f°(r 1 ,v 1 )=n- 1 J25°e(r 1 -q°,v 1 -w°), 
i=i 

where 5®(r,v) G Cq(G) (i.e., 5° is a continuously differentiable function in G with compact 
support), 5®(r,v) — >■ 8(r,v) and the support of 5® tends to {0,0} as e — > 0. The generalized 
function 5(r, v) is defined on the space T>(G) of test functions. Let \q® — q® \ > a for all % ^ j. 

1 ) Then, for every T > there exist a number et > and a family of functions 
D £ G CC 1 (fi), e > 0, such that there exists a solution f e (ri,vi,t) G C 1 (G x [—T,T]) of 
the Cauchy problem for the Boltzmann-Enskog equation with the initial data /°. This so- 
lution coincides with F lt£ (r 1 ,v 1 ,t) on (ri,Vx,t) G G x [—T,T] for all e < St- Here F 1>e is 
defined by the first formula of |3J) with D substituted by D £ . 

2) Consider an arbitrary time moment t G [— T,T] such that \qi(t,T) — qj(t, T)\ > a for 
all i 7^ j. Here qi(t,T), Wi(t,T), i = 1, . . . , iV ; and T are defined as before. Then, 

N 

lim/ £ (ri, v x ,t) = n^y^Jin - q i (t,T),v 1 - Wi(t,T)). 
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Proof. Consider the density function 



N 

O'lvJiv) i = 1 



where A = -u^ • • • > £am m at) £ G^, the sum is taken over all permutations ji, ■ ■ ■ , Jn of 
1, . . . , N. Due to the condition \q® — q® \ > a for all % ^ j and the conditions imposed upon 
S £ , for a given T, there exists So > such that D £ (A) = whenever |£° — < a for some 
2 7^ j, for all e < Eq. Hence, we can treat D £ as a function on Q 3 A. In this case, condition 
([5]) is also satisfied, since both sides of it are equal to zero. Note that by the same reason, 
condition ([5]) is also satisfied for the derivatives of D £ as well. Hence, D t£ = S_ t D £ also 
belongs to C 1 (f2) for all t and is continuously differentiable in t as well. 

Obviously, F 1 ° £ (r 1 ,v 1 ) = /°(ri, i^), where F® is defined by the first formula in ([8]) for 
t = and D substituted by D £ . 

For an arbitrary i G 1, we have 

1 N 

A,e(A) = S^D £ (A) = — U 5 e(^ A )-€M-t^)-<), 

where t, A),ui(— t, A), . . . t, A),un(— t, A) are the positions and momenta of N 

hard spheres whenever their initial positions and momenta were A. There exists E\ > (let 
£i < £o) such that the last formula can be expressed as 

1 N 

' O'lvJJv) * =1 

whenever e < ej. Here £*(r, t, T) — > 5(r, v) as £ — > for all z = 1, . . . , N and for all t such 
that \qi(t,T) — qj(t,T)\ > a, i ^ j. Intuitively: if e is sufficiently small and the points A 
with nonzero D(A) are concentrated near T, it is almost equivalent to shift A on —t along 
the dynamical flow or to shift r on t along the same flow. For a given T > 0, there exists 
St > (let St < S\) such that the support of 8l(r,v,t) lies in B a x M 3 for all i, t G [—T,T], 
and e < Et- Here B a = {(r, v) G G : |r| < a}. 

The functions J 7 ! and F 2 given by (jSJ) with D substituted by D £ have the forms 

N 

F h£ ( Xl ,t) = n~ l 6 l( r i - ii(t, r), Vl - Wi {t, r) ,t,r), 

F 2 , £ (x u x 2 ,t) = n- 2 (l - ^) _1 X)^( r i -Qi(t,r),V! - Wi (t,T),t,T) 



x5i(r 2 - qj {t,r),v 2 -w j (t,r),t,r). 



If s < e T , then 



„ . Fi £ (xi,t)F l £ (x 2 , t), if\r 2 - ri \>a, 
F 2 , £ {xi,x 2 ,t) = { . (11) 

0, otherwise 
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for all t G [— T,T]. 

The functions Fi and F 2 satisfy equation ([I]). But in view of (jTTj) . equation ([1]) is reduced 
to the Boltzmann-Enskog equation for Remind that F^ e (xi,0) = F® £ (xi) = f £ (xi). 
Hence, Fi ie (x\,t) is a solution of the Cauchy problem for the Boltzmann-Enskog equation 
with the initial data /° whenever e < e?- Fi jE G C'(Gx [—T,T]) in force of the contin- 
uous differentiability of D t£ (h) in t and A. The assertion of paragraph 1) of the theorem 
has been proved. The assertion of paragraph 2) immediately follows from the facts that 
f £ = F l £ and 5 l £ (r, v, t, T) — >• 5(r, v) as e — > for all i — 1, . . . , N and for all t such that 
\q i (t,T)-q j (t,r)\>a,i^j. □ 

So, if the density function in the Boltzmann-Enskog equation has the form of a sum 
of continuous delta-like distributions, then these distributions behave as "spreaded hard 
spheres": they collide with each other, move freely between the collisions, and gradually 
spread. The Boltzmann-Enskog equation ([1]) coincide with the first equation of the BBGKY 
hierarchy fflQj) . When the spatial spreading of the hard sphere centers exceeds the hard sphere 
diameter, equations ([T]) and ffTUl) cease to coincide. Note that this moment T increases as 
£ decreases. If e = and f e is a sum of delta-functions (j3j), then T = oo. This is the 
case of propositions 1 and 2. So, the approximate microscopic solutions turns to the exact 
microscopic solutions in the limit. 

The considered solutions f £ can be called the "approximate microscopic solutions" . Like 
the exact microscopic solutions, they are reversible, but only on the bounded time segments. 

So, the microscopic solutions of the Boltzmann-Enskog equation are not just a paraphrase 
of another equation fflOl) . Undoubtedly, the microscopic solutions reflect some fundamental 
properties of the Boltzmann-Enskog equation. They have the following meaning: they 
describe the limiting motion of the centers of delta-like distributions in the approximate 
microscopic solutions. 

Note also that equality ffTT]) is an analogue of the molecular chaos condition (or Stoss- 
zahlansatz) for the Boltzmann-Enskog equation. 

5 Reversible and irreversible classes of solutions of the 
Boltzmann-Enskog equation 

Now we are able to answer the question about the relation between the irreversible 
Boltzmann-Enskog equation ([1]) and its reversible solutions (jlj). For generality let us 
denote A the generalized functions in form (j4j) with arbitrary functions qi(t,T), Wi(t,T), 
i = l,. ..,N. 

Theorem 1.1) Let us treat the Boltzmann-Enskog equation (TJP as an equation for functions 
f G C 1 (G x [0,T)) ; T > (the classical solution of the Boltzmann-Enskog equation is 
defined in this class or its subclass 12^ [23f). In this case, the equation is irreversible, i.e., 
if some function f(ri,Vi,t) is a solution of this equation, then, in general, the function 
f(ri,Vx,t) = /(ri, —Vi, —t) is not a solution; 

2) Let us treat the Boltzmann-Enskog equation as an equation for generalized functions 
f G A on test functions from CC l (VL). In this case, the equation is reversible, i.e., if some 
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generalized function f(ri,vi,t;T) is a solution of this equation, then the generalized function 
f(ri,Vi, t; T) = /(r 1; —V\, —t\ V) is also a solution. 

3) The continuously differentiate solutions f E {fi,Vi,t) constructed in theorem 1 (sums 
of delta-like functions) are reversible on time segments [—T,T], where T > depends on 

£>0. 

Proof. Since = mn,-vu-t) = _ df^,-vi,-t) ^ WQ haye 

df(r 1} v u t) df{r u -v u -t) 2 f , , 
jr. = -Vi ^ a / \(v 21 , <r)|[/(ri, -v l9 -t)/(ri + aa,-v 2 ,-t) 

- f{n, -vx, -t)f(n - aa, -v 2 , -t)}dadv 2 . 

Replace f(r, —v, —t) = f(r, v, t) and change the variable a — > —a: 

dj{r u v u t) df(n,v l7 t) 2 f ~ ~ 

777 = ~vi ^ ha / {v 21 ,a)[f(r 1 ,v 1 ,t)f{r 1 + aa,v 2 ,t) 

at ar l J{V2l,T)>0 

-f(n, v[,t)f(n - aa, v' 2 , t)]dadv 2 . 

The Boltzmann-Enskog equation is obtained from this equation by swapping the velocities 
v\ )2 and v[ 2 in the collision integral. Hence, the function / satisfies the Boltzmann-Enskog 
equation if and only if 

f(ri,V!,t)f(ri - aa,v 2 ,t) = f(r u v[, t)f(r x - aa,v' 2 ,t), 

or, equivalently, 

/(ri, v u t)f(n - aa, v 2 ,t) = /(n, v[, t)f(n - aa, v' 2 ,t), (12) 
for all ri, Vi, a and t. 

Condition (1121) is, in general, not satisfied if / is a classical solution. Hence, if solutions 
from C l (G x [0,T)) are considered, the Boltzmann-Enskog equation is irreversible. 

Now let / 6 A. In this case, condition f|T2l is satisfied, since it is equivalent to (Q by 
definition. Indeed, both sides of f|T2|) are understood as functionals on CC 1 ^): 

J f(r 1 ,v 1 ,t;T)f(r 1 -aa,v 2 ,t;T)D(T)dT = J f^v'^t^f^-aa^t^D^dT. (13) 

This is exactly condition (Q. Hence, if solutions from A are considered, the Boltzmann- 
Enskog equation is reversible. 

Finally, if we consider the solutions / e (ri,t>x,t) constructed in theorem 1, then condition 
(|T2|) is also satisfied if t G [—T,T]. Truly, it coincides with ([9]) (with F 2 replaced by F 2)£ ), 
since f £ (xi,t)f £ (x 2 , t) = F 2>£ (xi, x 2 , t) for t G [-T, T\. □ 

So, condition (Q for test functions is crucial for the reversibility of the Boltzmann- 
Enskog equation for the class of solutions A. If we drop this condition, then, in general, the 
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microscopic solution flj]) does not satisfy the Boltzmann-Enskog equation. The right-hand 
side and left-hand side time derivatives of / are different, only the right-hand side time 
derivative is equal to the right-hand side of the Boltzmann-Enskog equation. Using a more 
sophisticated analysis we can prove the assertion of theorem 1 for this well. 

This result can be compared with the Kozlov's diffusion theorems for a collisionless 
continuous medium in a rectangular box [3]. The dynamics is described by the Liouville 
equation. If we consider a solution in the form of a delta-function, then the solution is peri- 
odic or almost periodic (since this case is reduced to the single-particle dynamics). But if we 
consider an integrable solution, the spatial density weakly converges to the uniform distri- 
bution, as time indefinitely increases. So, the asymptotic dynamical properties of solutions 
depend on the considered class of functions. 

Here we show that in the case of the Boltzmann-Enskog equation, not only the asymptotic 
properties of solutions, but also the reversibility or irreversibility property of them depend 
on the considered class of functions. 

6 Reconciliation of the kinetic equations with the mi- 
croscopic dynamics 

Of course, the above arguments cannot be regarded as a derivation of the kinetic Boltzmann- 
Enskog equation from the microdynamics. According to the traditional point of view, we 
must start with some initial iV-particle distribution function D (a randomization of the 
initial iV-particle phase point) and analyse the dynamics D t = S_pD of it. By doing this, 
we get equation f flO|) (as well as the other equations of the BBGKY hierarchy). Then we 
must obtain a kinetic equation using some limiting procedure. But, as we have stated in the 
introduction, until now we have not a fully satisfactory realization of this program. 

However, the microscopic solutions of the Boltzmann-Enskog equation and the inves- 
tigations of the reversibility or irreversibility property of it suggest us another way of the 
reconciliation of the kinetic equations with the microdynamics. 

Let us postulate not the microdynamics, but the Boltzmann-Enskog equation. That is, 
the dynamics of a system that we call "a system of hard spheres" is postulated to be described 
by the Boltzmann-Enskog equation ([[J with a continuous function /. Then, this equation 
is noticed to have also special (noncontinuous) solutions in form @ (or the continuous 
solutions constructed in theorem 1). The dynamics of ((&(£, T), Wi(t, T)), i = 1,...,N, in 
these solutions is interpreted as the microscopic dynamics. 

So, the establishment of the existence of microscopic solutions for a kinetic equation 
can be regarded as a derivation of the particle dynamics from it. Instead of the derivation 
of the kinetic equations from the microdynamics we suggest a kind of derivation of the 
microdynamics from the kinetic equations. 

According to this point of view, the "fundamental" laws of motion are not the laws of 
motion for individual particles (reductionism), but the laws of motion for the function /, 
which describes the system as a whole (holism). Laws of motions for the individual particles 
are approximations. This approximation is suitable if the density function has a special 
(delta-like) form (see, e.g., theorem 1). 
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Functional mechanics [TTJ [121 D2J EEH EE [16] also claims that the Newton equation (or 
Hamilton equations) is not fundamental even in the classical nonrelativistic world. 

Also, it should be noted that, in principle, the function / in the Boltzmann-Enskog 
equation, as well as other kinetic equations, can be interpreted without any reference to 
the notion of a particle. Let / be normalized not on the volume V, but on the total mass 
M = Nm (m is the mass of a particle): J f(r,v, t)drdv = M. Then the Boltzmann-Enskog 
equation has the same form, but the constant n = y in the collision integral ([2]) is replaced 
by — . f(r,v,t)ArAv has the meaning of the mass of the matter which is contained in the 
volume element Ar and moves with the velocity in the range Av. So, f(r, v, t) is the density 
function of the matter in the phase point (r,v) G T 3 x M. 3 . 

We have a kind of continuous medium. In contrast to the usual continuum mechan- 
ics, there is no definite velocity in a fixed point of physical space: in the same vol- 
ume element Ar different parts of the matter move with different velocities. The terms 
f(ri,Vx,t)f(ri ± acr,V2,t) in the collision integral are intensities of interaction of different 
parts of the matter. 

Initially, atomic theory of matter was an attempt to reduce the properties of the matter 
to the properties of its elementary blocks. But the problems like the reversibility paradox 
show that this program and the notion of an atom should be revised. 

Thus, we may still interpret the function / as the single-particle density function, or we 
may interpret it as the density function of a continuous matter. In both cases, we postulate 
the Boltzmann-Enskog equation for the system as a whole rather than the equations for the 
individual particles. 

Remark. A question concerns the notion of the derivation can arise. Mathematically, we can 
postulate a kinetic equations and speak about the derivation of the particle dynamics from 
it. But actually, kinetic equations are known to be obtained heuristically with the use of the 
notion of a particle and the laws of motion for particles. 

It is more preferable to speak not about a derivation, but about a reconciliation of the 
kinetic dynamics with the microdynamics. Usually, the reconciliation is understood as a 
derivation of the kinetic equations from the microdynamics. We claim that neither kinetic 
dynamics can be derived from the microdynamics, nor microdynamics can be derived from 
the kinetic dynamics in the "pure" sense. Every level of nature has its own laws, but one 
can speak about some agreement between them. 

We do not demand a kinetic equation to be derived from the microdynamics, by we say 
that a kinetic equation is in agreement with the microdynamics if it has microscopic solution 
or if it is an approximation of another kinetic equation which has microscopic solutions (see 
below) . 

7 Other kinetic equations 

7.1 Boltzmann— Enskog— Vlasov equation 

The obtained results can be generalized to some other kinetic equations. Bogolyubov [TJ [2] 
noticed that the kinetic equation for hard spheres with additional long-range pair inter- 
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action potential $(|r$ — rj\) also has microscopic solutions. Such system is described by 
the Boltzmann-Enskog equation with an additional "Vlasov integral" (we will call it the 
Boltzmann-Enskog-Vlasov equation) : 

777 = -Vi Vna I {V2i,cr)[f{ri,v 1 ,t)f{r 1 + aa,v 2 ,t) 

at ar l J(v2i,er)>0 

u + \f( a , n /" 9$ (l r l ~ r 2|) / i\ j 9f(ri,V!,t) 

- /(ri,vi,*)/(ri - aa,v 1 ,t)\dadv 2 -\ / p{r 2 ,t)dr 2 = , (14) 



m J dr\ dv\ 

where m > is a mass of a hard sphere, p{r 2 ,t) = f f(r 2 ,V2,t)dv2- Let us formulate the 
proposition analogous to proposition 2. 

Proposition 3. Consider the (generalized) function f given by (0) where the time depen- 
dences of qi(t,T) and Wi(t,T) are defined by initial values T and the Newtonian laws of 
motion: (ji(t) = Wi(t), and 

^) = ivWl« (15) 
W d qi {t) 1 ; 

if \Qi(t) — Qj(t)\ > a f or a tt 3 7^ h otherwise Wi(t) changes according to the law of elastic 
collision (with the replacement of v\ t 2 by Wij). This generalized function is defined as a 
functional on CC 1 (J1) in the same manner as in section 3 (with the correspondingly defined 
dynamic flow Sj: N ^). This function satisfies the Boltzmann-Enskog-Vlasov equation fLfy . 

A theorem analogous to theorem 1 can be formulated for this case as well. The formula- 
tion is exactly the same, except of the corresponding definition of the dynamic flow. 

Theorem 3. 1) Let us treat the Boltzmann-Enskog-Vlasov equation [Lfy as an equation 
for functions f G C 1 (M 3 x IR 3 x [0, T)), T > 0. In this case, the equation is irreversible, 
i.e., if some function f(ri,Vx,t) is a solution of this equation, then, in general, the function 
f{r\,V\,t) = f(ri,—Vi,—t) is not a solution; 

2) Let us treat the Boltzmann-Enskog equation as an equation for generalized functions 
f 6 A on test functions from CC 1 (fi). In this case, the equation is reversible, i.e., if some 
generalized function f(ri,vi,t;T) is a solution of this equation, then the generalized function 
f(ri,Vi,t;T) = f(ri,—vi,—t;T) is also a solution. 

3) The continuously differentiate solutions f E (ri,Vi,t) constructed in theorem 1 (sums 
of delta-like functions) are reversible on time segments [—T,T], where T > depends on 

£>0. 

The proof is analogous to the proof of theorem 2 (the Vlasov term is reversible in both 
cases, only the Boltzmann-Enskog collision integral is irreversible in the case of a continuous 

/)• 

Thus, we can speak about a derivation of the classical Newton equation (or Hamilton 
equations) from the Boltzmann-Enskog-Vlasov equation, since it has the corresponding 
microscopic solutions. 
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7.2 Boltzmann equation for hard spheres 

Let us consider the Boltzmann equation for hard spheres. It has the same form ([1]) as the 
Boltzmann-Enskog equation, but the collision integral has the form 

St / = na 2 / (v 2 i, o-)[f(ri,v[,t)f(ri,v' 2 ,t) _ f(r 1 ,v 1 ,t)f(r 1 ,v 1 ,t)}dadv 2 . 

J(v2i,er)>0 

The Boltzmann equation has not microscopic solutions. It can be obtained from the 
Boltzmann-Enskog equation by the Boltzmann-Grad limit a — > 0, n — > oo (or A" — > oo, 
since n = y and V = const), na 2 = const (also one can prove the week convergence of 
the corresponding solutions [2U 125]). The Boltzmann-Enskog equation is considered to be 
more exact than the Boltzmann one [19]. The Boltzmann equation is an approximation of 
the Boltzmann-Enskog one in the case of a large number of hard spheres and a negligible 
diameter of them. 

It is worthwhile to note that the Bogolyubov's derivation of the Boltzmann equation from 
the microdynamics leads not to the Boltzmann equation itself, but to some modification of 
it, which turns to the Boltzmann-Enskog equation in the case of a hard sphere gas [5]. 

So, the Boltzmann equation has not microscopic solutions, but it can be represented as 
an approximation of an equation which has such solutions. In this sense one can speak about 
a reconciliation of the reversible microdynamics of hard spheres and irreversible dynamics 
of the Boltzmann equation for a hard sphere gas. 

7.3 Boltzmann equation in the general form 

Consider the Boltzmann equation with a more general form of the collision integral 

7TT = -Vi k + na / B(v 21 ,a)[f{r l ,v 1 ,t)f(r 1 ,v 2: t) 

ot dn J R3xS 2 

- f{n, i>i, t)f(r 1 ,v 1 ,t)]dadv 2 . 

This case describes the dynamics of the particles with a pair interaction potential $. The 
function (collision kernel) B is determined by $. 

This equation also has not microscopic solutions. However, as in the considered particular 
case of hard spheres, it would be interesting to obtain this equation as a limit of another 
kinetic equation which has microscopic solutions. Let us sketch this construction. 

Let $ be an interaction potential which corresponds to the given B. Let TV particles 
with positions qi(t) and velocities Wi(t) at the moment t, i — 1, . . . , N, interact with the pair 
potential $( ^'~ 9 ^ ), a > 0. Assume that $(r) -» as r — > oo. 

Consider the Boltzmann-Grad limit — > oo, a — > 0, Na 2 — > const. Since a — > and 
$(r) — > as r — > oo, the interaction time tends to zero. Hence, the interactions can be 
approximately described as instantaneous. The dynamics of particles with instantaneous 
interactions can be described by a Boltzmann-Enskog-like equation, which has microscopic 
solutions. 

A detailed construction of this approximation will be a subject of a further work. 
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8 Conclusions 



We have formulated the existences of microscopic solutions of the kinetic Boltzmann-Enskog 
equation for hard spheres and the Boltzmann-Enskog- Vlasov equation for hard spheres 
with additional long-range interaction potential as rigorous propositions and have pointed 
out the importance of additional condition fl5]). We have analyzed the relation between 
the irreversibility of these equations and the reversibility of the microscopic solutions of 
them. The solutions of these kinetic equations are time-irreversible if the considered solutions 
belong to the class of continuous functions. They are reversible if the considered solutions 
have the form of sums of delta-functions (microscopic solutions have this form). So, the 
reversibility or irreversibility property of the dynamics depends on the considered class of 
functions. 

Also we have constructed the continuous approximate microscopic solutions. They are 
reversible on bounded time intervals. The approximate microscopic solutions have the form 
of sums of delta-like distributions. They behave as "spreaded hard spheres", the limiting 
motion of the centers of these delta-like distributions is described by the exact microscopic 
solutions. 

We have used these results to propose a new way of reconciliation of the kinetic equations 
with the microdynamics. In contrast to the traditional paradigm, where the laws of motion 
for individual particles are regarded as fundamental laws (reductionism), we have postulated 
the kinetic equations, which describes the dynamics of the system as a whole (holism). The 
establishment of the existence of microscopic solutions for a kinetic equation is a kind of 
derivation of the microdynamics. So, one can speak about a derivation of the microdynamics 
from the kinetic dynamics. 

Finally, we have discussed the Boltzmann equation. It has not microscopic solutions, but 
it can be regarded as an approximation of the Boltzmann-Enskog (or a Boltzmann-Enskog- 
like) equation in the case of a large number of particles and a negligible interaction radius 
of them. In this sense one can speak about a reconciliation of the microdynamics with the 
Boltzmann equation as well. 
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